Dephasing of spin and charge interference in helical Luttinger liquids 
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We consider a four-terminal Aharonov-Bohm interference setup formed out of two edges of a 
quantum spin Hall insulator, supporting helical Luttinger liquids (HLLs). We show that the tem- 
perature and bias dependence of the interference oscillations are linked to the amount of spin flips in 
tunneling between two HLLs which is a unique signature of a HLL. We predict that spin dephasing 
depends on the electron-electron (e-e) interaction but differently from the charge dephasing due to 
distinct dominant tunneling excitations. In contrast, in a spinful Luttinger liquid with SU(2) invari- 
ance, uncharged spin excitations can carry spin current without dephasing in spite of the presence 
of e-e interactions. 
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Quantum spin-Hall (QSH) insulators can support edge 
states that are topologically protected against disorder, 
and are present in the absence of a magnetic field PHI] • 
These edge states have a special helical structure, in that 
their direction of propagation is associated with a given 
value of spin polarization. Evidence of the existence of 
such helical edge states has been recently found in HgTe 
based quantum well structures [3 [6]. They in princi- 
ple enable direct control of spin currents via electronic 
means, which has generated interest in their spintronics 
applications [7J, [S] - 

However, the fact that spin can be controlled electron- 
ically implies that it must also couple to the electron- 
electron (e-e) interaction. This is in contrast to SU(2) 
symmetric systems in which propagating spin excitations 
can be uncharged, and therefore remain unaffected by 
Coulomb forces [5]. The e-e interaction is the domi- 
nant cause for dephasing of electronic coherence in low- 
dimensional mesoscopic systems, at temperatures low 
enough [10], and in ID systems, this dephasing is asso- 
ciated with fractionalization of excitations [TT]. In QSH 
systems, one would therefore expect that dephasing orig- 
inating from the Coulomb interaction would be visible 
also in the spin coherence, and lead to fractionalization of 
spin [12] . On the one hand, this will limit the spin coher- 
ence time which may be harmful for some applications, 
but on the other hand, it can serve as a characteristic 
signature of the QSH state. 

One possible way to probe the coherence in a low- 
dimensional mesoscopic system is to observe the decay of 
interference oscillations |13] . A simple geometry in which 
this can be done is an Aharonov-Bohm (AB) interfer- 
ometer, where the phase difference between two paths of 
propagation is controlled with a magnetic flux. Although 
such flux would not ordinarily couple to spin currents, 
the structure of helical edge states makes this possible 
[7]. Alternatively, one can observe the interference os- 
cillations as a function of gate or bias voltages, or spin 
density differences, applied to the system — these will be 
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FIG. 1. (Color online) Interferometer supporting helical edge 
states, where the direction of propagation is correlated with 
the spin (t / i). The four-terminal setup is connected to 
noninteracting leads, biased at potentials fti/2, l/r- Tunnel- 
ing between the edges occurs at points x = and x = d 
(shaded). The loop between the tunneling points is threaded 
by a magnetic flux <1>. 



present also in non-helical edge states. 

Here, we analyze dephasing of charge and spin interfer- 
ence in an AB interferometer composed of tunnel-coupled 
interacting helical edge states (see Fig. [T]) . We find that 
dephasing in charge tunneling currents depends on the 
degree of spin-flips in tunneling. Moreover, we estab- 
lish a direct relationship between charge and spin cur- 
rents in helical four-probe setups, which translates the 
results on charge dephasing to apply to spin. Here, heli- 
cal edge states turn out to differ significantly from ordi- 
nary spinful edge states, which in the presence of SU(2) 
spin symmetry can support interference effects in spin 
current without dephasing, even when Coulomb interac- 
tion is present. 

Model. An interacting helical edge state can be mod- 
eled as a helical Luttinger liquid (HLL) [3]. Each HLL 
has the same degrees of freedom as a spinless LL (SLLL) 
with right (+) and left (— ) movers associated with a given 
spin state (|+) =|t)>| — ) =14-}); so s P m is redundant. Mi- 
croscopically, the spin states |f) and |i) refer to Kramers 
partners which are either electron spin or total angular 
momentum eigenstates, depending on the material. 
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The bosonized Hamiltonian for a HLL reads [3J 

/oo 
Ax[g{x)- 2 : {d x df : + : {d x 4>f ■] (1) 
-oo 

where the standard boson fields $(x), 4>{x) satisfy 
[<j){x),'&{x')] — (i/2)sgn(x — x') and are associated with 
annihilation operators for left- and right-going electrons, 
ip a {x) oc e MkFX+^{x)]+iV^4>(x) ; a = ±j and g ( x j ig 

an e-e interaction parameter. Moreover, vf, kp are the 
Fermi velocity and Fermi wave vector. For a pair of edge 
states, the spin-direction mapping is reversed on the sec- 
ond edge. We assume the edge states are contacted to 
nonintcracting Fermi leads a distance I away from the 
central region. That is, there is a repulsive interaction 
{g{x) = g < 1) in the central region (— I < x < I), but the 
leads (\x\ > I) are noninteracting (g — 1). Here, we as- 
sume that the relevant lead modes coupled to the system 
are also described by Eq. 0. We will use dimensionless 
units in which h — e = kg = I. 

The tunneling between the upper and lower edges is 
described by the Hamiltonian 



(It) via the Kubo approach (cf. [131 [2U]). The tunneling 
current operator can be identified as It = /+ + /_, where 



h t = l ]T M^V-i,, 

R=0,da/3=± 



! (#)tV> 2 , /3 (i?)+h.c, (2) 



where j = 1,2 denote the upper and lower edges and L is 
the total length of an edge (including the leads). Under 
time-reversal symmetry broken only by the magnetic flux 
$ in the interferometer loop, the tunneling elements read 
t ap (R) = e-^W^/*H> ap (R), if ap = -aW^y, 
where $o = h/e is the magnetic flux quantum. In the 
presence of inversion symmetry, spin is conserved [2] and, 
as a consequence, the spin-flipping terms should vanish, 

t' ++ = t' ~ 0. However, local gates [Hj or strain can 

induce spin-mixing in the tunneling amplitudes, allowing 
for non-zero t' aa . 

We compute currents in the presence of tunneling, 
treating Hj- as a perturbation. This is expected to be 
valid in HLLs, as at low energies the tunneling stays ir- 
relevant (in the renormalization group sense) for 1/2 < 
g < 2 [T3 US]. Bias voltages in terminals are taken 
into account by assuming that the incoming states have 
thermal populations described by chemical potentials 
/ii/2,± = A*i/2,L/fl an d a temperature T. These potentials 
can be gauged into the tunneling Hamiltonian [IT] [T5] : 
t a p(R) >-» t aP (R)e l ^ la -^ 2l3)t in the interaction picture. 
That the system is contacted to Fermi leads implies that 
despite any fractionalization to counterpropagating plas- 
mons, all charge injected to the + (— ) channel finally 
enters the right (left) lead [TS]. Moreover, we assume 
that I is large compared to the interferometer size d and 
length scales h/T, h/V given by the bias and tempera- 
ture, so that we can neglect any finite-? effects [18] on the 
tunneling dynamics. 

Within this approach, the total tunneling current from 
the upper edge to the lower edge is found by computing 



iL 22 2J t a/3 (R)ip ly . 

R=0,dP=± 



t (R)ifo,p(R) + h.c (3) 



describe tunneling into the + and — channels. 

Charge tunneling current. The leading contribution to 
the charge tunneling current is given by single-particle 
processes, as described by H T , in the range 1/2 < 
g < v3 [HI [TB]. In the loop geometry, there are 
three distinct ones: direct tunneling current through each 
contact separately (\t a /3(x)\ 2 ), an interference contribu- 
tion involving both contacts (t a< - a (0)t a ,-a(d)*), and 
an interference contribution with spin flips in tunneling 
{t a ,a(fi)t a ,a(d)*), the latter two in general coupling to 
the AB phase. In the leading order in tunneling, particle 
conservation forbids mixed processes (e.g. i ++ t + _). 

For later convenience, we first consider the current 
tunneling into the a = ± channels separately, (I a ) — 
jdirect _|_ jab _|_ /AB,sf T nese are gj ven \yy the expressions 



7 direct = ^|^( J R)| 2 Re^ 1 (0,/ii, a -M2^) (4a) 

R=0,dl3=± 

l£ B = 2 Re[t a ^ a (Q)t a ^ a (dYe- 2lak - d ] (4b) 
x ReZ 2 (d, n ha - (i 2 - a ) , 
7 AB,sf =2Re[t aa (0)t aa (d)*Z 1>a (d,^ a -ti 2 , a )} (4c) 

where Z\_ + = Z\ _ = Z\, with Z\(x,ui) — 
tL 2 JZ o dte^ t ([(4^ 2 , + )(x,t) 1 (4 i+ ^ h+ )(0,0)}) Q and 
Z 2 (x, u>) = iL 2 e 2lk ^ dt e™* { [(^,+^2,- ) (x, t) , 
(ip2 _ Vi.+)(0, 0)]) Q . These equilibrium correlators of the 
clean system can be found via bosonization techniques 
[5], and the time integrals can be evaluated analytically. 

We concentrate on the tunneling charge current in a 
situation where a bias V is applied between the upper 
and lower edges, fii. a = V/2, (i 2 . a — —V/2. At high 
temperatures (ttT/u = z 3> l/d), we observe that the AB 
oscillations with flux <I> experience exponential dephasing 
(see Fig. [2]), which for V <C T gives 



Jab « Z 2 



4 2 ^2nn 2 F v 2 F z(azrsin^ h+1/2)zd 



(5a) 

4 2 ~<2irn F v 2 F {az)^ ivd/u iyzd 

JAB,sf oc Zi ~ e ' Ve (5b) 

(1 + 2j)u z 



Here 7 = \{g + g x ) - \, 



VF/g is the plasmon 



velocity, and rip = L/(2itvf) the density of states at the 
Fermi energy, and a the short-distance cutoff length. 

The finite exponent <i"fzd occurs due to charge frac- 
tionalization [T3J. At low temperatures, V <C T <C u/d, 
the exponential dependence crosses over to a power law 
in T, as shown in Fig. [2] and Z\ and Z 2 coincide. 
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FIG. 2. (Color online) Temperature dependence of 
the amplitude of the dimensionless conductance corre- 
sponding to interference oscillation amplitudes, G = 
(27rn|) _1 Re &Z 1/2 /&V \v=o- The spin-flip G A B, s f (solid), 
non-spin- flip Gab (dashed), and the high-temperature scal- 
ing from Eq. |5| (dotted) are shown. 




FIG. 3. (Color online) Bias dependence of the amplitude 
of the interference component of the current, for T = 0. 
The spin-flip (solid), non-spin- flip (dashed) contributions are 
shown. Apart from the low-bias power law, the noninteract- 
ing (7 = 0, g = 1) and interacting (7 = 0.05, g w 0.64) cases 
differ mainly in the oscillation period nu/d = irvp/gd. 



Note that for Tab , the oscillations dephase exponen- 
tially even in the absence of interactions (7 = 0). The 
difference arises because spin-conserving tunneling in a 
HLL fixes the direction of propagation, and so the pair of 
trajectories contributing to interference has an unavoid- 
able dynamical phase difference of 2kd which also leads 
to a dependence on kp [21] [Jab oc cos(2fcpd + ipo)] [22] . 
similarly to what happens in chiral Luttinger liquids |20j . 

Moreover, the interference contributions in the current 
oscillate not only with the flux, but also with the bias V, 
as shown in Fig. [3j The single characteristic frequency 
iru/d is different from interference in a SFLL, where two 
characteristic frequencies exist [HI [23] because of spin- 
charge separation. 

Relation between charge and spin currents. To un- 
derstand the dephasing in the spin current, it is useful 
to note first that in a HLL the spin current is closely 
associated with the charge current. The changes in 
charge (c) and spin (s) currents due to tunneling are 
5I C oc Sp + — 5p- , and 6I S cx 5p+ + Sp_ in terms of 
changes in the densities p. Note that this spin current is 



defined with respect to Kramers partners. When evalu- 
ated in the noninteracting right (left) lead, in which the 
+ and — modes are independent and Sp_ = (5p + = 0) 
because of ballistic transport, we find 5I^ R = 8I] R and 
5ll L = —5I] L on the upper edge 1 (positive direction 
of current is from left to right in Fig. [T]). On the lower 
edge 2, the sign of the spin current is flipped. In a four- 
terminal setup, this relation becomes more transparent 
if one considers an XYZ decomposition [TS] of the cur- 
rents, extended to account for possible non-conservation 
of spin: I x = \ [I 1L + I 2L + I 1R + I 2R ] indicates the 
total current flowing from the left (L) to the right (R), 
I Y = \ [I lL -I 1R -I 2L +I 2R ] the total current from top to 
bottom, I z = \[I 1L + I 1R -I 2L -I 2R ] the current flowing 
between the diagonals, and I s = -I 1L - I 2L + I 1R + I 2R 
the non-conserving "source" current flowing out of the 
system. In this representation we find: 



sir 

si? 




(6) 



since the dc charge current is conserved (7f = 0). The 
result applies to any helical four-terminal setup. Note 
that the spin current flows perpendicular to the charge 
current, which is a signature of the spin-Hall effect giving 
rise to the helical state. As a result, the spin current can 
be accessed by a measurement of the transverse charge 
current. 

Spin tunneling current. From Eq. ^ it immediately 
follows that also spin currents suffer from e-e interaction 
induced dephasing. For instance, consider the spin tun- 
neling current (Y) that can be generated by applying a 
bias V s between left and right (X), pi + = /12+ = V s /2, 
III- = P2- = —Vs/2- In this configuration, there is no 
(Y)-charge current [24l . similarly as in the (Y)-bias con- 
figuration where there was no spin (Y) -current. Here, 
the spin current is given by the difference between parti- 
cle currents tunneling into the + and — channels, and 
it can be directly evaluated using Eq. Q: (It,s) = 
(-I+ + L) = If ilcct + 7 S AB . To leading order in H T , 
the spin (Y)-tunneling current has the same form as the 
charge (Y)-tunneling current in the Y-biasing configu- 
ration (see Eq. (5a)), except that in this order in Ht, 
spin- flipping tunneling cannot contribute [25j . 

Depending on the interaction parameter g and the 
amount of spin-flips in tunneling, the interference contri- 
bution in spin current can at high temperatures be dom- 
inated by two-particle tunneling processes (see Fig. |4| , as 
described by the effective Hamiltonian (cf. |15j). 

i4 2) = v p i/; 1+ 4_xp 2+ 4_ (7) 
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(a) 



Vp,sf 




(b) 



V p , s { 



FIG. 4. (Color online) Dominant (neutral) compound tunnel- 
ing contributions to spin interference (Y)-current in (X)-bias 
configuration, (a) Illustrations of compound processes. The 
process v p conserves spin, whereas process v PiS f flips one spin, 
(b) The dynamical phase differences of interfering paths for 
these compound terms with arrows describing the propaga- 
tion directions of electrons [filled dots in (a)]. 



The spin-density fluctuation assisted backscattering pro- 
cess Vp^f, which is present only when spin flips are al- 
lowed, has a slower dephasing than single particle pro- 
cesses in the whole range 1/2 < g < 1: IP' s f' AB oc 
V a cos(V s d/u) cos(2kpd + ip ) exp(— 2gzd). The spin- 
conserving process v p on the other hand dominates 
the single-particle one only for g < l/y/3: If ,AB oc 
sin(2V s d/u) cos{Akpd + ipo) exp(—4gzd). Although these 
processes do not transport charge between the two edges 
nor couple to the flux $, the oscillation with the bias V s 
and hp remains. 

Note that the above implies that in a four-terminal 
HLL setup, charge and spin tunneling currents have dif- 
ferent dephasing exponents, as the interference contri- 
bution is dominantly carried by different types of exci- 
tations: Neutral electron-hole excitations [Eq. Q] for 
g < 1 have different correlation lengths compared to 
single-particle excitations. Consider, for instance, the 

single-particle dephasing exponent [1 + (^-§ l)]2zd 

from Eq. ( 5a I 



compared to the v p process, Agzd = 
[1 — (1 — g)]Azd, where 2zd and Azd arise from dynami- 
cal phase differences [cf. Fig. Qb)]. Interestingly, for v p 
the repulsive interaction reduces the dephasing from the 
nonintcracting value. 

The above should be contrasted to what occurs in usual 
SFLLs. There, cotunneling of electrons can create an un- 
charged spin excitation (spinon), which does not couple 
to interactions. The lowest-order tunneling process cou- 
pling only to the spin sector is the two-electron process 



H. 



(2s) 



X 1 



T ll>l,a,-a +h.C 



(8) 



where opposite spins tunnel to opposite edges. A calcula- 
tion along the same lines as above yields the interference 
component in the spin tunneling current, 

7 s , osc = 47rr4K| 2 V;'[(V;') 2 + (2vrT) 2 ] cos(V^d/v F ) , (9) 



where V' s = /iif — fJ-ii — /i2t 4" ^21 ■ Exponential dephasing 
is indeed absent, as expected under SU(2) spin symme- 
try. Generating such spin currents requires a "spin bias" 
Vg, which in our model is essentially equivalent to a dif- 
ference in spin densities between the two edges. In SFLLs 
one possibility for inducing this is to contact the SFLL 
to a system in which a spin imbalance is externally main- 
tained 26J, or to couple it to a HLL 27\ which may also 
allow measuring the spin currents via charge currents. 

In summary, we showed how electron-electron inter- 
action results in dephasing of interference oscillations in 
ID helical liquids, both in charge and spin tunneling cur- 
rents, with respective exponents that can differ. More- 
over, we pointed out how the close coupling of the spin 
current to the charge current in a helical liquid can re- 
sult in a qualitatively different behavior from spinful Lut- 
tinger liquids. Such effects provide a clear signature of 
the helicity of the transport, and understanding them 
may be valuable for applications. 

We thank B. Trauzettel and T. Ojanen for useful dis- 
cussions, and acknowledge financial support from the 
Emmy-Noether program of the DFG. 
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Current operators 

In this work we are interested only in the time-averaged 
spin or charge currents entering the noninteracting leads. 
To compute them conveniently, we need to identify the 
operators corresponding to the tunneling currents, taking 
into account the plasmon reflections at the edges of the 
leads. 

The identification can be done using a similar approach 
as in Ref. [19] . From the Heisenberg equation of motion 
for the field </>, under a Hamiltonian H — Hq + V[4>, 
where the effective tunneling term V has support only in 
a finite region, one finds an exact result for the change in 
the charge current operator SI C — n~ 1 / 2 VFd x ((j) — <f>\v=o) 
in the Heisenberg picture: 



SI, 



/CO 
dt 'Ax ' 
-°° a=± 



n:| — — -D a (xt, x't') 

29 



1-9 



29 

]a[X> 2J^\S6(x' 



D- a {xt,x t )jj a (x t ) , 
SV 



+ a- 



2^ \<50(x') 5d{x')J ' 



(10) 



(11) 



where D± are solutions to the plasmon wave equation 
d 2 D — d x (v F g(x) d x D) = with initial conditions cor- 
responding to right (left) moving (5-pulses starting at 
x = x': D±(x,t',x',f) = 5(x - x'), d t D± (x, f, x' , t') = 
=pu6'(x — x'). Reflections at the lead edges where g(x) 
changes split the initial pulse to a train of pulses escap- 
ing to the leads. Solving the wave equation, one can 
verify that the fraction of an initially right-going pulse 
entering the right (left) lead in the long-time limit is 
w R+ = lim^oo J °° d(t - f) D + (x,t,x' ,f) = (l + .g)/2 



1+gl+g 1-gl-g 

~Yg 2^ 2^ 

1 + 9 1 - 9 _ 1 - 9 1 + 9 

~Yg 2^ Hg 2^ 



(12) 
(13) 



Eq. (10 1 leads to the conclusion that (i) all current in- 
jected in the + (— ) channel enters in the right (left) lead, 
in the time average, and, (ii) I a = J dx' j a (x') is the op- 
erator corresponding to the total current injected into 
channel a. 



Several other points can also be directly read off ( 10 1 



since the fluctuation of spin density in a helical liquid 
is proportional to the same operator 8p + — <5p_ as the 
charge current, it can be seen that a spin injected to a 
helical liquid fractionalizes into counterpropagating plas- 
mons carrying the fractional spin (g ± l)/2g (cf. 12 ). 
A similar calculation for the charge density fluctuation 
implies that an injected electron fractionalizes to plas- 
mons carrying the charge (1 ± g)/2, as usual 28J. The 
equation also guarantees that similarly as charge, a spin 
injected to the + (— ) channel will, in the long-time limit, 
be transmitted in entirety to the right (left) lead. This 
implies that spin fractionalization cannot be detected by 
simple time-averaging spin current measurements (such 
as those proposed e.g. in [H]), but instead one needs to 
have access to time scales of gl/vp characteristic of the 
plasmon transport. 



Scaling dimensions and dephasing 

The limiting behavior of dephasing of interference ef- 
fects at low or high temperatures can be found from an 
analysis of the scaling dimensions of the terms in the 
effective tunneling Hamiltonian. 

Given a bosonized operator M = 

ve tn^(4,+§) e iq^(<j>-8) ( Qn a sing l e edge ) ; w ithin 

our model one finds the correlation function 



(M(x,*)M(0,0)t) oc \v\ 2 F(x,t) 2A +(F(x,-t)*) 2A - , 

(14) 

A ( \ [(l-9)"+(l + 9)g] 2 a i , ne ,s 
A+(n,g) = =A_(g,n), (15) 



where F(x,t) = iazj sinh[z(x — ut — ia)], z = irT/u, and 
a is the short-distance cutoff length. The scaling dimen- 
sion of operator M is A = A + + A_. Now, if M repre- 
sents an effective tunneling term in the Hamiltonian, a 
similar Kubo calculation as done below indicates that at 
low temperatures its contribution to spin/charge conduc- 



tance scales as G 



but at high temperatures 



Td/u ^> 1 interference effects dephase as G oc e ~ 47rTdA <i ; 
where the dephasing exponent A^ = min[A + ,A_] can 
differ from A/ 2. 
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Comparing between effective tunneiing processes 
allows us to identify those dominating at high temper- 
atures, and results for the scaling dimension A indicate 
the magnitude of the tunneling element, and justify the 
use of perturbation theory. 

For 1/2 < g < 2 we have A > 1 for all time-reversal 
symmetric processes, so that in view of renormalization 
group flow, the situation is perturbatively stable, and 
tunneling scales to zero at low energies [TJl US]. Con- 
cerning dephasing of interference effects, of the tunnel- 
ing processes transporting charge, we find that at high 
temperatures the single-particle process dominates (for 
.9>l/2), 

A d -(.g + .g- 1 )/4-l/2. (16) 

If spin flips are not allowed, the dominant process is still 
a single-particle one 

A d = ( ff + ff - 1 )/4, (17) 

and the situation stays the same in the whole range 1/2 < 
g < 2. 

For spin transport, the situation is considerably dif- 
ferent: suppose first that spin is conserved. Then, for 
l/\/3 < g < >/3 single-particle tunneling dominates, 

ipi+4-, A d = ( 5 + . 9 - 1 )/4, (18) 

but for 1/2 < g < a two-particle cotunneling pro- 

cess has the lowest exponent, 

V>i+4-V>2+</4_ , \i = g- (19) 

In this case, it occurs that A^ = A/2. 

If spin flips are allowed, one finds that the process with 
the smallest A^ in the whole range l/2<g<lisin fact 



a two-particle one, 

^ 1+ V4+V>2+V4- , A d = g/2 . (20) 

Note that (i) this process is not important for the charge 
tunneling current, since it does not transport charge, and 
(ii) it is higher order in tunneling and has a larger scal- 
ing dimension A = 1 + g than the single particle process. 
Note that it can also be written in the explicitly time re- 
versal symmetric form ipi + ip\_ [^4+^2+ — V4-V'2-] + h.c. 
resembling spin density fluctuation assisted backscatter- 
ing. 

From the above discussion, we conclude that the lead- 
ing results for charge tunneling (Y) current are obtained 
with first-order perturbation theory, but for spin current 
also second-order contributions, or alternatively, the ef- 
fective 2-particle tunneling, needs to be analyzed. 



Kubo correlators 

The correlation functions Z\, Z 2 appearing in the 
Kubo calculation are 

/oo 
dte^K+OM^+OM), (21) 
-oo 

4 + (0,0)^i, + (0,0)]_) , 

/oo 
dte iw *([^ )+ (x,t)^ J -(aj,t), 
-OO 

(22) 

4_(0,Q)Vi J+ (0,0)]_) . 

These correlators can be evaluated via standard 
bosonization techniques for spinless Luttinger liquids [5] 



Z\(x,u) — UnpVpa 2 
Z2(x,uj) — 2inpV F a~ 



— oo 
oo 



dte iut Im 
dte ilJt Im 



2 7 +2 



shxh(z(x — ut + idj) i Vsinh(z(a; + ut — id)) 



2 7 +l 



sinh(z(a; — ut + idj) > Vsinh(z(a; + ut — id)) 



(: 



2 7 1 



27+1" 



(23a) 
(23b) 



Here, a is the short-distance cutoff, 7 = j(g + g^ 1 ) — |, 
u = vp/g is the plasmon velocity, z — nT/u the inverse 
thermal length, and np — L/(2itvf) the density of states 
at the Fermi energy. 



Integrals of the above type can be evaluated in closed 
form by making use of the binomial series 

( 1 \ _ 2 a -a(sgnHeO r ( a + U ) -(a+2wXsgnR«flg 

(24) 

and its Fourier transform. In the limit a = Im£ — > + 
one finds 
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dte ~ T7 x ( 25 ) 

> \smh(w — t) J 

2°tt 2 iru/2 

r(a) e 

T(l + f - f )r(l - f - f )[cosh(vrw) - cos(Tra)] ' 

which has simple poles at uj — ia + 2m, n = 0, 1, 2, 



We now rewrite (23) in the form f(uj) — /(— w)*, where 
/ is a convolution that can be transformed via contour 
integration to a Matsubara sum, which can be evaluated. 
As a result, we find 



Zi{x, U ) - 427+W ^ (aZ)47 * ( fg^e-rd + 2 7 - ^) 2 F l( 2 7 , l + »y-p-, (26) 
u \1 (2 + 27) 2uz 2uz 2wz 



-4(7+l)xa 

i^r^ 1 + 27 + 2^ 2 + ^ 1 + 27 + ^ 2 + » e 



RcZ ^ <"> = u r ( i + 2 7) Im [ g r ^ + 27 - j*^* 1 + 1 + 2 ^ 2W 1 - 2W 6 )] ' 

(27) 



where 2-Pi(a, c, z) = 2-Fi(a, c, z)/T(c) is the regular- 
ized hypergeometric function. For z> l/rf and V <C T 
we obtain Eqs. (5a) and (5b) of the main text. 



Compound tunneling between HLLs 

The spin transport due to the compound tunneling 
Hamiltonian 

i4 2) = v p ^ 1+ 4_ij 2+ 4_ (28) 

can be handled similarly as above. For Ht — 
v p ipi+ip\_'i])2+' l p2- + h.c. we find I + = —I- = 

2ivpil>i + t/jl_ip2+*l>2- + h.c. on the basis of Eq. (111. From 
this we immediately see that in the X-bias configuration 
the charge (Y) current vanishes, and the spin (Y) current 
is obtained by a straightforward calculation: 



Similarly, starting from 



(It,s) = 2 



E I' <" 

R=0,d 



v p (R)\ 2 ReZ 3 (0,-2V s ) 



(29) 



+ 4Re[v p {0)v p (d)*e 4 * kFd ] Re Z 3 (d, -2V.) , 
Z 3 (x,u) = (2iray 2 Z 2 (x,uj)\ 1 ^ g _ 1 /2 ■ (30) 

Apart from the prefactor, the result is identical to Z2 



4 2) = v- 



(31) 



we get I 1+ = = 2iv PiS fip 1+ ipl_[ij}l + ip2+ - 

ipt-tpz-} + h.c. = iM\ — iM\, and similarly on edge 2. 
The relevant correlation functions are of the form 

(M 1 (r,t)M 1 t (0,0)) o oc |«p )Bf | V < <*« 1 +-' 41 ->*e 2i ** r (32) 
x F(r,ty(F(r,-ty) 9 [F(r,t) 2 + (F(r,-t)*) 2 ] , 

which implies that the high-temperature dephasing expo- 
nent is Ad — g/2, being different from A/2 = g/2+1. We 
remark that the contribution in third order in tunneling 
vanishes identically due to particle conservation. 

Note that in the X-bias configuration, we have 81% = 
27+, 51% = 0, and 51% = for both of the above pro- 
cesses. 



Description of SFLL 

In a spinful Luttinger liquid [5], the direction of spin 
and motion is not coupled and the fermion operator has 
to be defined separately for spin and direction of motion 
(at each edge): 



with a different exponent, in agreement with Eq. (14) 



(33) 
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It is advantageous to change to charge (p) and spin 
0) variables T(i) = (l/V^Op ± <f>a) and # t( j.) = 
(l/>/2)(0p±i?a), explicitly 



a iak F x pi-v/f (0p+Qi9 p +s(0< T +Qi?„)) 



(34) 



with s = ± for a =^, \.. In these new variables, the in- 
teracting system generally splits into a sum of spin and 
charge parts: i? S FLL = H p + H CT +g 1 dx cos(v / 87r^ tT ). 
These Hamiltonians are characterized by interaction pa- 
rameters g v and velocities v v , v = p, a. At the low-energy 
fixed point, g\ flows to zero for repulsive SU(2) invari- 
ant interactions and therefore the Hamiltonian becomes 
H = H SLL = H p + H a with 



and 



H a = — I dx 



9 P {dxfipf + — {d x -dp) 2 
9p 



9a (d x (j)a) 2 + — {dx^crf 

9<y 



, (35) 



(36) 



where g p = g, u p = u, and g a = 1 [9]. The spin sector 
therefore becomes effectively non-interacting. 

The dephasing for spin currents can therefore be absent 
in second order in the tunneling which we describe by the 
effective Hamiltonian 

cr=t4- Q 

(37) 

e.g. consider the term 

i> 1+t 4 +t rP 2+l 4 u ~ e V^Wi*+**-fc*-*»«>. (38) 

The correlation function of terms like this one will not 
show fractionalization because g a = 1 and the charge 
fields have disappeared from the bosonized expression. 
The coupling constant v a in the effective Hamiltonian 



Eq. ( 37 ) is second order in the bare tunneling amplitude 



and can depend (at low bias voltages) on temperature 
in a power-law fashion (see section below) but will not 
depend on the interferometer length d as long as Td/u 
1. 



Renormalization of tunneling operators 

For completeness, we now summarize how compound 
tunneling terms can be derived via the real-space per- 
turbative renormalization group |29j . which has proved 
useful in studies of Luttinger liquids [30] ■ 

We first note that bosonization allows evaluation of 
correlation functions involving Fermi operators in closed 



form. In general, for a set of bosonic operators <&j 
well-known expansion applies: 



(39) 



where T indicates (contour-)time ordering. Given 
the bosonic correlation functions, ($±(i)$±(0))o = 
^- g ln{-izL/ smh[z{t - ia)}} for $± = 2~ 1 / 2 (</> =p $/#), 
one finds a short-time expansion (TA(t)B(t')Q)o oc \t — 
t'\- a (TA{t)B(t)Q) a as t' -> t for groups A and B and 
Q of Fermi operators. Since the scaling at t 1 — > t is 
independent of Q, in view of Eq. (39), this can be un- 



derstood as an operator product expansion A(t)B(t') = 
\t-t'\- a (AB)(t) ast' ^t. 

The perturbative renormalization group now proceeds 
by integrating out all short-time divergences at times 
a + da > \t — t'\ > a appearing in the second order 
of perturbation expansions. The generated compound 
terms are absorbed in the first order of the expansion via 
a change in the Hamiltonian, H h-> H + dH , which after 
rescaling the cutoff a H> a + da leads to the RG flow 
equations. 

In this approach, we find the flow equations corre- 
sponding to the effective Hamiltonian £f c g = Hq + Ht + 



H, 



(2) 



= [l-(g- x +g)/2]t aP , 



dt a p 
din a 

= [2 - 2g]v p - cat+^+Oig- 1 - g) , 

d In a 

^ = [1 - 9]v P . si - c'at ++ t*_ ^{g- 1 - 1) 
d In a 



(40) 
(41) 
(42) 



where the first terms appear, as usual, from the intrinsic 
cutoff dependence of the operators, ■ e** J )o oc a A , and 
the latter from the short-time divergences. The factors 
c and c' are numerical constants. Note that the effec- 
tive Hamiltonian is written in terms of Fermi operators. 
Moreover, the compound terms above are generated only 
for repulsive interactions, g < 1 — for g > 1 the pertur- 
bation expansion does not have the corresponding short- 
time divergences. 

The flow equation can be integrated starting from the 
bare values a = a 0l t a p(a a ) = t a p, v p (a ) = v p ^{a ) = 0, 
up to the length scale at which the operator product 
expansion breaks down. In our case it is the size of the 
interferometer d or the thermal length 1/z, whichever is 
smaller. This yields the scaling of the prefactors with the 
cutoff (assuming 1/z <C d): 

t^^iaoT^+^n-h^O), (43) 
v p oc (aor) 111 ^ 29 - 1 ^ -1 ^- 2 ^- 1 ^^^)^ + (0) , (44) 
v p , s{ oc (a T) min (^^ 1 +f- 2 )T- 1 t ++ (0)r + (0) . (45) 

A similar treatment gives the scaling for the process de- 
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scribed by Eq. (37) in a SFLL 



« (aoT) min ( 1 '( g ~ 1+g ^ 2 ~ 1 )T~ 1 t + f !+ f(0)t+, + ,(Q) , 



(46) 



The result for v p coincides with that obtained in Ref. [TS] 
for g < (v5 — 1)/2 at which the scaling 2g — 1 of the com- 



pound process starts to dominate the scaling g^ 1 + g — 2 
of two single-particle events. The dephasing exponents 
for the above processes, on the other hand, follow from 
the arguments in the previous sections, and in general 
they are independent of the scaling dimensions of the 
prefactor. 



